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ABSTRACT 

The total cosmic-ray luminosity of the Galaxy is an important constraint on models of 
cosmic-ray generation. The diffuse high energy 7-ray and radio-synchrotron emissions of 
the Milky Way are used to derive this luminosity. The result is almost two orders of 
magnitude larger than the standard estimate, based on the observed isotopic abundances 
of cosmic ray nuclides. We discuss the plausible interpretation of this discrepancy and 
the possible origin of such a relatively large luminosity. 
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Almost a century after they were discovered, our understanding of cosmic rays is still 
very limited. Their production mechanisms, composition and energy spectrum continue 
to be debatable. In this letter we discuss the total cosmic ray (CR) luminosity of our 
Galaxy, a crucial constraint on models of galactic CR generation. 

The CR nuclei have a power-law spectral flux dF/dE oc E"'''' with a series of break- 
point energies: Pi^2,3 ~ 2.7, 3.0, and 2.5 in the intervals 10 GeV < E < Eknee ~ 3 x 10^ 

GeV, Eknee < E < E^nkle ~ 3 X 10^ GcV, and Eankle < E < 3 X W GeV. Below Ek„ee, 

protons constitute ~ 96% of the CRs at fixed energy per nucleon. Their fiux and number 
density above Ep ~ 10 GeV is (see, for instance, Wiebel-Sooth and Biermann 1998 and 
references therein): 
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The spectral indices of heavier nuclei are compatible within errors with that of protons 
and, given the dominant abundance of the latter, we need not distinguish here between 
CR protons and the ensemble of nuclear CRs. It is generally believed that the bulk of 
the CR nuclei with energy below the knee are Galactic in origin, and that their main 
production mechanism is acceleration by supernova shocks (see, for instance, Ginzburg 
and Syrovatskn 1964; Longair 1981; Berezinskn et al. 1990; Gaiser 1990). 

If the CRs are chiefly Galactic in origin, their accelerators must compensate for the 
escape of CRs from the Galaxy, in order to sustain the observed Galactic CR intensity: 
it is known from meteorite records that the CR flux has been steady for the past few 
giga-years (Longair 1981). The Milky Way's luminosity in CRs must therefore satisfy: 

LcR^Lp = -/J-E^dEdV, (2) 

C J Tconf dE 

where Tconf (E) is the mean conflnement time in the Galaxy of CRs of energy E. 



The standard estimate of Lcr runs along the following lines. The mean column density 
X traversed by CRs before they reach the Earth can be extracted from the observed ratios 
of primary to secondary CRs. The result is X « 6.9 [E(GeV)/(20 Z)]~°-^ g cm"^ (Swordy 
et al. 1990). With use of X = /pdx ~ pc Tconf one can extract the product of Tconf (E) 



and a path-averaged density p. Assume the locally measured values of X and dFp/dE to 
be representative of the Galactic values dominating the integral in Eq. (Q), to obtain: 

L„.l^/pdv/iE|^dE. (3) 

Assume the path-averaged p to be close to the average density p of neutral and ionized 
gas in the Galaxy, so that / p dV is the total mass of Galactic gas, estimated from X-ray, 
optical and radio observations (Longair 1981) to be Mgas ~ 4.8 x 10^ Mq. The integration 
over energy is not unduly sensitive to its lower limit and converges rapidly above the knee. 
The final result (Drury et al. 1989) is: 

LcR ~ 1.5 X 10^^ erg s"^ . (4) 

Earlier estimates (e.g. Berezinskii et al. 1990 and references therein) of Lcr, which used 
the "Leaky Box" model of CR confinement, led to somewhat smaller luminosities. 

In spite of the cursory character of the above luminosity estimate, Eq. (P is consistent 
with the assumption that CRs are dominantly accelerated by the turbulent magnetic 
fields of supernova (SN) remnants, generated by the expansion of the debris from the 
SN explosion into the interstellar medium. For an estimated mean Galactic rate of one 
supernova every ~ 50 years (van den Bergh and Tammann 1991) and an average kinetic 
energy (E^) ^ 10^^ erg of the debris, this explanation requires an e ~ 20% efficiency in 
the conversion of kinetic energy into CRs. 

The model of CR generation by SNe is incomplete or problematic in several respects 
(e.g. Plaga et al. 1999). Supernova-generated shocks are not sufficiently lasting and 
energetic to produce CRs with energies well above Eknee (e.g. Lagage and Cesarsky 
1983). The space distribution of SNe is too concentrated in the Galactic disk and bulge 
to give a proper description of the relative isotopic abundances of CRs, in particular 
^°Be/^Be (Strong and Moskalenko 1988), of the directional distribution of the diffuse 
7-ray background radiation (Strong and Mattox, 1996) and of the high energy 7-rays 
produced by CR interactions in the interstellar medium (Strong and Moskalenko 1998). 
The CR luminosity of SN remnants is severely constrained by TeV 7-ray observations and 
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results in a CR-generation efficiency e between 1 and 5% (Allen et al. 1999), somewhat 
short of the required e ~ 20%. 



In this letter we present an alternative estimate of the luminosity of CR nuclei, based on 
the Galactic CR-electron luminosity, which we infer from observations of 7-ray production 
(Hunter et al. 1997; Sreekumar et al. 1998) and synchrotron emission (e.g. Chen et al. 
1996) by CR electrons. 

The EGRET detector on the Compton GRO sateUite has mapped the intensity and 
spectral index of the "diffuse" 7-ray background (GBR) above = 30 McV, at latitudes 
above the Galactic disk and bulge. The observed spectrum, dF/dE^ oc E^^"', has an index 
(3^ ~ 2.10 ± 0.03 that is independent of direction. The intensity is also roughly isotropic, 
thus the claim of a dominantly extragalactic origin of the GBR. 

We have recently shown (Dar and De Riijula, 2000) that the GBR intensity is signif- 
icantly correlated with the angle away from the galactic centre and that it is dominated 
at high latitudes by inverse Compton scattering (ICS) of CR electrons from the cosmic 
microwave background radiation (CBR) and from starlight, obviating the recourse to un- 
specified extragalactic sources (the importance of ICS has also been stressed in , Strong 
and Moskalenko 1998; Strong et al. 2000). Earlier evidence for a large galactic contribu- 
tion to the GBR at large latitudes had been found by Chen et al. (1996), who discovered 
a strong correlation between the observed EGRET GBR 7-ray intensity and the galactic 
radio continuum emission at 408 MHz, which is dominated by synchrotron radiation from 
the very same CR electrons that produce ~ 100 MeV 7-rays by ICS from galactic stellar 
hght. 

Our model of the origin of the GBR is based on the assumption that the average 
Galactic CR-electron spectrum has the same energy dependence as the locally observed 
one (for a recent compilation of experimental results see Wiebel-Sooth and Biermann 
1998). This spectrum is well fit, from ~ 10 GeV to ~ 2 TeV, by: 
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Starlight and CMB photons, upscattered by electrons with the spectral index of Eq. (j^), 
have an energy dependence with index [3^ = (/5e + l)/2 = 2.10 ±0.05, in perfect agreement 
with the observed direction-independent index of the GBR (Dar et al. 1999, Dar and De 
Riijula, 2000). 

The relation between the spectral indices of CR protons and CR electrons in Eqs. (|l]) 
and d^) can also be understood in very simple terms. The confinement, residence or 
accumulation time of nuclear CRs in the Galaxy is rconf(E) oc X oc E^°-^. The result must 
be of the same form for electrons, since relativistic particles of the same charge behave 
in the same way in a magnetic maze. The source spectrum of nuclear CRs, dFp/dE, is 
related to the observed spectrum of CR nuclei by dFp/dE oc rconf(E) dFp/dE, so that the 
index of dFp is (3^ = jSp — 0.5 ~ 2.2. If the mechanism accelerating CR hadrons and CR 
electrons is the same e.g., first-order acceleration by a moving magnetic field (Fermi 1949, 
1954) dFp/dE oc dF^/dE, and the source spectral index for electrons is also (3^ ~ 2.2. 
Electrons, unlike nuclei, are significantly affected by ICS and synchrotron cooling -whose 
characteristic time is TcooI oc m^/E- so that, at sufficiently high energy, cooling takes over 
the accumulation-time effect in modulating the electron spectrum. The result (Dar et al. 
1999, Dar and De Riijula, 2000) is /3e = + 1 ~ 3.2, in agreement with Eq. (§. 

In our study of the GBR we adopted, for the spatial distribution of the CR electron 
flux in the Galaxy, a model with a gaussian scale height he above the Galactic plane, and 
a scale radius pe in directions perpendicular to the Galactic axis. By adjusting hg ~ 20 
kpc and pe ~ 35 kpc, we reproduced the observed intensity and angular dependence of the 
GBR. We shall assume the nuclear CRs to be distributed as the CR electrons, with the 
ratio of fluxes fixed at its locally observed value. A scale height of CR nuclei hcR = he ~ 
20 kpc is larger than conventionally assumed, but is not excluded by data on relative CR 
abundances. For the most elaborate models (Strong and Moskalenko 1998) a "Leaky-Box" 
scale height of 20 kpc is only some 1.3 standard deviations below the central value of the 
most precise observations (Connell 1998) and is perfectly compatible with the average of 
all previous and somewhat less precise results, compiled in Lukasiak et al. (1994). Since 
our CR distribution is much more extensive than the visible part of the Galaxy, we refer 
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to it as the "cosmic-ray halo" . 



The EGRET GBR data to which we fit the properties of a CR electron halo are 
gathered by masking the galactic plane at latitudes |b| < 10°, as well as the galactic 
centre at |b| < 30° for longitudes |1| < 40°. The volume of the CR halo is so much larger 
than that of the Galaxy within EGRET's mask, that it is a very good approximation 
-in computing the Galaxy's total CR luminosity- to use our halo model throughout the 
entire Galaxy (within the mask the model accounts for ~ 1/2 of the observed diffuse 7 
radiation) . 

The successful relation between the spectral indices of the GBR and the CR electrons 
followed from the assumption that the production rate of CR electrons is equal to their 
cooling rate, which we estimate as follows. The starlight-photon density in the CR halo 
may be obtained by approximating our galaxy's starlight as that produced by a source at 
its centre with the galactic luminosity = 2.3 x 10-'^'^ ~ 5.5 10^^ eV s~^ (Pritchet and 
van den Bergh 1999): n* ~ L^/(47rce*r^), where e^, ~ 1 eV is the average photon energy. 
For a gaussian CR halo, the mean n^, is given by: 

(n,) P. - In (l^) ^ 0.035 cm"^ , (6) 

47rce^p; u Vl — uy 

with u^ = 1 — hg/pg. The mean energy density of starlight in the CR halo is much smaller 

than that of the GBR (~ 0.24 eV cm~^). If the local magnetic field-energy is in equipar- 

tition with the CR energy density, B^/Stt ~ 1 eV cm~^, the mean total electromagnetic 

energy density in the halo is ^ 1.27 eV cm~^. For the electron energy range of interest 

the Thomson limit of the 67 cross section (a^ ^ 0.65 x 10~^^ cm^) is accurate, even for 

ICS on starlight, and the mean cooling rate, Rc, of CR electrons by ICS and synchrotron 

radiation is (Dar and De Riijula 2000): 

The luminosity of our galaxy in high energy electrons of energy above E, in equilibrium 



1 4p^(7„cE 
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with their cooling rate by ICS and synchrotron radiation, is: 



By substituting the flux of Eq. (|^), we obtain: 



Lc(> E) ^ 1.13 X 10 
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We liave assumed that the ratio of the CR nuclear and electron fluxes is universal 
throughout the Galaxy. Thus, to estimate the Galaxy's CR luminosity, it suffices to 
scale the electron luminosity Lg by the local ratio R of CR and electron fluxes. For the 
total fluxes, this ratio is R ~ 80. This result is uncertain, since it is dominated by CR- 
energies of 0{1) GeV, a domain in which the fluxes are affected by local magnetic and 
solar-wind effects. An independent estimate of R can be obtained as follows. The CR 
electron spectrum sharply steepens to an index jSe — 3.2 at E ~ 5 GeV: that must be 
the energy at which ICS and synchrotron radiation, for which TcooI oc E^^, take over the 
effect of CR accumulation, for which Tconf oc E~°'^ (the effects of local magnetic fields, the 
solar wind. Coulomb scattering, ionization losses and bremsstrahlung are only relevant at 
even lower energies). Thus, to within a factor of 0{2), the observed proton to electron 
flux ratio at E = 5 GeV must be the ratio of their source fluxes. We have argued that 
the source fluxes have the same spectral index. Thus, their ratio at fixed energy is also 
their energy-integrated ratio. This gives R ~ 60, in rough agreement with the previous 
estimate. 



Multiply Le in Eq. (g) by R = 60, to obtain: 



Lp(> E) ~ 6.8 X 10 
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which is almost two orders of magnitude larger than the estimate of Eq. 



If, as we argued, Tcooi = Tconf at E ~ 5 GeV, we can use Eq. (^ to obtain Tconf = 250 My 
at E = 1 GeV. This is an order of magnitude larger than the values of Xconf obtained from 
the analysis of the relative abundances of unstable to stable CRs: ^^Be/Be (Lukasiak 
et al. 1994; Connell 1998), ^^A\/^'^A\ (Lukasiak et al. 1994b; Simpson and Connell 
1998) and '^^Cl/Cl (Connell et al. 1998). The discrepancy is even larger, since these 
data are for lower energies, of (9(250) MeV per nucleon. This alterity can be easily 
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understood (Plaga 1998). The confinement time is extracted from tlie isotopic ratios 
using a Leaky Box model, wlierein tlie magnetic field of the Galaxy is confined to a region 
of dimensions similar to those of the visible part of the Galaxy. But, if the dense and 
luminous component of the Galaxy is embedded -as we surmise- in a much larger and 
less dense magnetized halo, the stable CRs may spend much of their travel time in the 
halo, while the unstable ones must have much shorter trajectories (the lifetimes of ^°Be, 
^^Al and ^^Cl are a mere 1.6, 0.87, and 0.30 My, respectively). This is also the reason 
why our estimate of the CR luminosity, Eq. (|10D , is not truly contradictory to the much 
smaller conventional estimate of Eq. : the gas density, volume and grammage used to 
derive Eq. all refer to CRs confined to a region close to the visible Galaxy. 

What CR acceleration mechanism could give rise to the large luminosity of Eq. (jlOD? 
The bulk of the high energy CRs may be accelerated by relativistic jets emitted in the birth 
of neutron stars and stellar black holes in supernova explosions (Dar and Plaga 1999). If 
the mean sky velocity of neutron stars (Lyne and Lorimer,1964), (vns) — 450 ± 90 km s~^, 
is due to an imbalance in this relativistic jet ejection, Ejet > Mns ^ns c ^ 4 x 10^^ erg, for 
= 1.4 Mq. If the kinetic energy of the jets is efficiently converted to CR energy, 
and for the estimated rate (van den Bergh and Tammann 1991) of Type II, lb and Ic 
supernovae, Rsn ~ 1/50 per year, then: 

LcR ~ 2 Ejet Rsn ~ 5.1 x 10^^ ^-i ^ (H) 

in reasonable agreement with Eq. (plQl). 

Relativistic jets are emitted by active galactic nuclei and by galactic microquasars. 
These jets are observed to consist of "plasmoids" whose cross section, after an initial 
period of transverse expansion at the speed of sound in a relativistic plasma (c/v^), 
remains surprisingly constant until the jet sweeps enough material to stop and disperse 
as a blob (Rodriguez and Meribel 1999). If the plasmoids of the jets allegedly responsible 
for the peculiar velocities of neutron stars have a Lorentz factor 7 = E/Mc^ of O{10^), 
they are good candidate sources of 7-ray bursts (Dar and Plaga 1999). If the transverse 
size of these plasmoids of C(O.l) pc, the column density necessary to stop them is of the 
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same order as the one transverse to the galactic disk. Thus, the jets may reach the halo 
of the Galaxy before they stop, seeding it with a CR population and a magnetic field, 
and giving consistency to our overall picture. 
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